clc;
clear;
A=[16 2 3 13;
    5 11 10 8;
    9 7 6 12;
    4 14 15 1];
p1=poly(A),%求解特征多项式
norm(p1-[1,-34,-80,2720,0])
p2=poly1(A),%求解特征多项式,利用的是书本的p109的poly1函数